Nonlinear management of the angular momentum of soliton clusters 
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We demonstrate an original approach to acquire nonlinear control over the angular momentum 
of a cluster of solitary waves. Our model, derived from a general description of nonlinear energy 
propagation in dispersive media, shows that the cluster angular momentum can be adjusted by 
acting on the global energy input into the system. The phenomenon is experimentally verified in 
liquid crystals by observing power-dependent rotation of a two-soliton cluster. 
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Introduction. - Angular momentum (AM) is a fun- 
damental quantity, the importance of which has been 
highlighted in almost all areas of physical sciences. Evi- 
dence of its role is found at the inception of the universe: 
although the distribution of galaxies, stars and planets 
is still a puzzle in astrophysics, it appears that an ini- 
tial angular momentum in the early universe prevented 
cluster-sized clouds from collapsing into a series of black 
holes, i.e. with no planets to support life [H, 0). In both 
classical and quantum mechanics, angular momentum is 
at the basis of rotational dynamics; hence, AM governs 
the behavior of important processes such those arising 
from fluid motion (e.g. initiation of cyclones Q and fluc- 
tuations in the length of a day Q), statistical complexes 
of rotating molecules [5(, quantum optics [6| and quan- 
tized particle ensembles. In the latter area, research has 
been mostly conducted in two major directions of inves- 
tigation: 

i) the revolution of trapped particles by their interaction 
with classical fields carrying angular momentum, foster- 
ing nanotechnological applications such as rotors or more 
complex machines powered by light [7[; 

ii) the preparation of energy packets in well defined AM 
states, with implications to both fundamental physics [8] 
and quantum information systems 

However, despite the importance of AM in physics and 
the vast literature on the subject, neither methods to 
nonlinearly control AM have been proposed, nor thor- 
ough studies have been carried out on the effect of non- 
linearity on the angular momentum of a specific system. 
Conversely, great attention has been devoted in the past 
few years to solitons and solitary waves. Such waves 
are ubiquitous and rely on the balance between wave- 
packet dispersion (spreading) and nonlinearity Fol- 
lowing the pionee ring numerical experiments by Fermi, 
Pasta and Ulam the universal concept of soliton 

has acquired importance in several physical sciences, in- 
cluding biology [l2| , hydrodynamic s lldl , plasma physics 

gravitation [ItJ 



condensates (BEC) — and are the subject of vigorous 
theoretical studies [l8J. Recently, the interactions be- 
tween two-dimensional deterministic soliton clusters have 
attracted interest, encompassing a wealth of dynamics 
ranging from wave filamentation (l9| , to spinning [20| ] to 
spiraling [2l|, [2J 



231 ] - Up to date, however, studies fo- 



14j . optics 



[131 ] . ultracold atoms 

and beyond [l(| . By virtue of their robustness, solitons 
have potentials in applications — from optical telecom- 
munications to atomic interferometry in Bose-Einstein 



cused on either specific nonlinear models or particular 
input waveforms 3, 21, 22, 23|, or on ID stochastic dy- 
namics [30, ]3l|; none of them discussed the role of the 
excitation on the dynamics of a two-dimensional multi- 
soliton ensemble. 

In this Letter, by employing a rather universal model 
for the theory and a nonlocal dielectric for the experi- 
ments, we investigate the behavior of a cluster of (2+l)D 
optical solitons, demonstrating a nonlinear approach to 
controlling its angular momentum. For the sake of sim- 
plicity but with no prejudice on its general validity, we 
develop the analysis in the simplest case of a two-soliton 
cluster, starting from first principles and employing the 
language of symmetries [241] . We demonstrate that the 
angular momentum of a soliton cluster exhibits a linear 
dependence on the nonlinear response; hence it can be 
precisely managed by varying the input power. The AM 
can therefore be evaluated from the global revolution of the 
cluster which, for a fixed propagation distance, evolves 
linearly with excitation. We check the theoretical re- 
sults by a series of experiments in nematic liquid crystals 
(NLC), a nonlinear nonlocal medium known to support 
stable (2+l)D solitons in voltage-tunable configurations. 

SI, Sill. 

Theory. We build a general model of nonlinear 

wave propagation, stemming from conservation laws and 
variational symmetries (as there is a one-to-one corre- 
spondence between them [24| ) and considering an iso- 
lated medium with translational and rotational invari- 
ance. The medium exhibits an optical response nonlinear 
with the wave intensity; hence it enforces the conserva- 
tion laws of momentum (M), angular momentum (A), 
Hamiltonian (H) and energy flux ( W) , respectively. Such 
a nonlinear model is derived by defining a suitable action 
integral T — J d 2 rdzC, the Lagrangian of which supports 
the variational symmetries originated by the basis of Lie 
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generators {vi = d/dx,v 2 — d/dy,v 3 — <9/<9z,v 4 = 
xd/dy — yd/dx,v 5 = iif) a d/dip a + c.c.} (we adopt Ein- 
stein's summation over repeated indices), being r = [x, y] 
and z dimensionless coordinates and ip a the dimension- 
less wave- function of the a— th wave-packet (a £ [1, 2] in 
this work). A general form for the Lagrangian L is as 
follows: 



' dz 



c.c 



x vcvv>" + d<n 2 e a/3 ® \^\ 2 (i) 



with a,0 e [1,2], V = [d/dx,d/dy], 9 Q/3 the Hermi- 
tian tensor and <g> a convolution operator defined by 
f®g — JJ d 2 r'f(r'— r)g(r'). The general character of the 
Lagrangian |T]) can be proven from the Euler-Lagrange 
equations of motion, which read: 







(2) 



The linear portion of Eq. ([2]) is Schrodinger like and de- 
scribes wave-packet propagation in the presence of dis- 
persion [2^]; the nonlinear portion {Q a p ® \4>^\ 2 ) can 
model both local — i.e. for 6 Q( g(r) = c a p5(r) — [HI 
and nonlocal fLU I27I |29| responses. By construction, 
the Lagrangian fl} admits the symmetries generated by 
Vj (i — [1,5]); hence Eq. ([2]) possesses the following in- 
tegrals of motion: 

M = J J <Fr{iip*Vip a + c.c.) (3) 

U =\jj d 2 r(\r\ 2 e«p ® l^l 2 - V^Vf) (4) 
A = J J d 2 r[r x (iip* a Vi> a + c.c.)] (5) 
W = JJ d 2 r\t(; a \ 2 (6) 

The conserved quantities ©-(JB]) allow to generalize the 
Ehrenfest theorem which, with reference to Eq. |(3J), takes 
the form: 



0_ 

dz 
0_ 

dz 



(i/j a \r\ip a ) = (ip a \p\ipa) 

(V'alPlV'a) = (V'alVGalV'Qt) 



(7) 



being Q a = Q a /3®\ip \ 2 , P = — «V the momentum opera- 
tor and (/alga) = // d 2 rf*g a the inner product defining 
the space metric. The set (0 is a generalization of models 
either derived in one dimension [301 ] or obtained through 
linear expansions 21, 23, 32]. Out of the whole spectrum 



of solutions to Eq. ([2|), we are interested in the evolution 
of a nonlinear cluster, the solitary "particles" (solitons) 
of which are found as invariant solutions of Eq. Q (writ- 
ten for a = (3) with respect to the global symmetry group 



generated by v = v 3 + v 4 - 
method of characteristics 
of each solitary wave: 

<>p a (r;z) = <j>{ 



,vi +y Q v 2 +^ v 5- Then the 
| yields the functional form 



r a \)exp(i/iz) 



(8) 



being [x a ,y a ] = r a = (i/j a \v\i/) a ) / JJ d 2 r\ip a \ 2 the soliton 
"center of mass". The term "particle" is justified as long 
as the solitons do not collide: the nonintegrable nature of 
Eq. ([2]), in fact, does not guarantee that after collision(s) 
the soliton ensemble evolves iso-spectrally (33|, i.e. with 
a constant number of wavepackets. However, if solitons 
simply overlap, the effects of the interaction terms (8 a /3<8> 
\4>/3\ 2 for a ^ (3) can be treated adiabatically, i.e. by 
assuming that the latter do not affect the functional form 
of the solitons but only their phase /1 or their center of 
mass r Q . Such hypothesis allows us to cast Eq. ((7|) in 
a potential form. By substituting Eq. JS]) into 0, after 
some straightforward algebra, we obtain: 



Q , ^ 







(9) 



being m — JJ d 2 T\%jj a \ 2 the soliton "mass", q = r x r 2 = 
[q x ,q v ], V, = [d/dq x ,d/dq y ], U(q) = -[U 12 (q) +U 21 (q)] 
and: 



Uat 



d 2 r>(|r' + q|) 2 e Q (|r'|), a ^ (10) 



Equation ^ describes the classical motion of a mass 
point (of mass m) subject to the nonlinear potential de- 
fined by (fT0| . It is worth to remark that we do not 
rely on a linear expansion; hence, Eq. ^ holds true 
also when the excitation varies. The Lagrangian density 
of ^ is C = qq/2 — U(q) and admits the variational 
symmetries generated by the basis {vi = d/dz,V2 = 
q x d/dq y — q y d/dq x }, the latter originating the conserva- 
tion laws of energy £ and angular momentum L, respec- 
tively: 



c qq 

t = 777 



L = TTiq X q 



(11) 
(12) 



Noticeably, equation (TT2"]) states that: 

i) the angular momentum is conserved, hence it is con- 
trolled by the system input; 

ii) the angular momentum depends linearly on the soli- 
ton "mass" 777, hence, on the nonlinear potential U. 

As a consequence, a change in the soliton global power 
JJ d 2 r\tp\ 2 = JJ d 2 r(|i/ji| 2 + |i/>2| 2 ) results in a linear 
variation of AM, provided the adiabatic condition [i.e. 
Eqs. ([8])-([TT])] holds. To measure the AM observable, 
hereby we suggest an original approach exploiting soli- 
ton spiraling. Such dynamics occurs when the initial 
momentum q(0) = fi(0) — £2(0) balances the attrac- 
tive force provided by the bound potential t/(q) [2l[, re- 
sulting in a rigid rotation of the (two) soliton ensemble 
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FIG. 1: (Color online). Sketch of the planar cell with nematic 
liquid crystals: (a) front view and (b) top view with an indi- 
cation of the molecular orientation (optic axis) in the plane 



with constant separation d and invariant angular veloc- 
ity lj — C/JVl — dS/dz, being S the angle spanned by 
q, M. the "momentum of inertia" of the classical system 
® [Hi . Following a straightforward integration of 
C/M. = dS/dz, a measure of soliton cluster AM is given 
by the revolution angle 8 which, for a given z, varies lin- 
early with the angular momentum C. As a result, accord- 
ing that the adiabatic condition holds, we expect that the 
revolution angle 5 of the soliton ensemble evolves linearly 
with the soliton input power P; n — JJ d 2 r\ip\ 2 . 

Experiments. In the experiments we used a pla- 
nar glass cell with a 100/im thick layer of planarly an- 
chored E7 liquid crystals (Fig. [I}. This configuration 
is similar to that previously employed for the investiga- 
tion of accessible solitons in highly nonlocal media 27j . 
but with a pre-orientation ^ rj 30 with respect to z 
in the (y, z) plane to make a voltage bias unnecessary 
[3~il | . We carried out the experiments with a near infrared 
(A = 1.064/^m) source and high resolution CCD cameras 
for imaging both the soliton output profiles and their 
propagation in (y,z). At the input, two extraordinary- 
wave solitons were excited with opposite momenta along 
x, compensating for the walk-off in (y,z). Particular at- 
tention was paid to the tilt in (x, z) in order to achieve 
spiraling with (invariant) separation d between the two 
"particles" . 

The experimental results are summarized in Fig. [2j show- 
ing images of light propagating in (y, z) (Fig. [5]a-d) and 
output intensity profiles in (x, y) (Fig. [5] e-h) for increas- 
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FIG. 2: (Color online). Summary of the experimental results: 
(a-d) evolution in the plane (y, z) and (e-h) output intensity 
profiles of the two-soliton cluster for increasing input powers 
Pin G [2.1, 2.7, 3.3, 3.9]mW. The borders of the NLC cell are 
indicated with a solid line in (e-h). 



ing excitation P m . In agreement with our theoretical 
prediction, as the power increased from 2.1 (Fig. [H^e) 
to 3.9mW (Fig. Eft-h), the cluster AM changed as well, 
as witnessed by the rigid rotation of ~ 180 degrees in 
the output plane (Fig. [2^-h). Remarkably enough, each 
soliton profile remained nearly unmodified (Fig. [3K) and 
the rotation angle S evolved linearly with power (Fig. 
[3})), demonstrating the nonlinear control over the overall 
AM angular momentum and in excellent agreement with 
the model. Owing to the giant reorientational nonlin- 
earity of NLC, the resulting sensitivity A<5/APj n = ir/2 
(rad)/mW is quite substantial. 

Conclusions. - Stemming from first principles and 
without specific assumptions on the dielectric response, 
we theoretically disclosed and experimentally demon- 
strated an original approach to gain nonlinear control 
over the angular momentum of a cluster of (two) spa- 
tial solitons. Such nonlinear management of the soliton- 
interaction potential is a remarkable example of all- 
optical control over light-induced guided-wave intercon- 
nect. Due to the general character of the model (and the 
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FIG. 3: (Color online), (a-b) Measured output spot-size w 
of the two solitons normalized to »o = w(2.1mW) and (c) 
revolution angle 5 versus input power Px n . 



experimental results) we derived, our findings are going 
to affect other areas where solitons are actively investi- 
gated, including plasma physics and Bose-Einstein con- 
densates, including novel applications in nonlinear optics. 
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